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INTRODUCTION 


The parameters which govern the mechanical behavior of materials also 
affect their stress wave propagation characteristics. In particular, stress 
waves are the interrogating energy in ultrasonic testing (UT) and acoustic 
emission (AE) . 

A major technique in UT is "through-transmission" ultrasonics. A 
mechanical disturbance is introduced into the material by a transducer coupled 
to one face of the structure and the resulting signal produced by a second 
transducer which is coupled to the directly opposite face of the structure 
is analyzed. Wave velocities, material attenuations, and structural thick- 
ness can be measured using this technique [1]. Also, cracks can be detected 
in this manner [2]. And, it has been shown that the inital through trans- 
mission longitudinal attenuation can be correlated with the compression 
fatigue life [3] and the flexural fatigue life [4] of graphite fiber 
composites . 

Another technique is "pulse echo" ultrasonics. Here the transmitting 
and receiving transducers are mounted on the same face of the structure. 

The output signal is due to reflections from the other face or a dis- 
continuity inside the material. The transmitting and receiving sensors may 
be separate or the same transducer. 

An ultrasonic parameter in which the transmitting and receiving 
transducers are coupled to the same face of the structure is called the 
"Stress Wave Factor" by Vary et al. [5,6]. An input pulse having a broad- 
band frequency spectrum is applied to the transmitting transducer and the 
number of output oscillations exceeding a certain threshold amplitude is 
defined as the stress wave factor (SWF). The SWF has been correlated with 
mechanical properties of graphite fiber epoxy composites [5-7], For 
example, Williams, and hamper t [7] showed that as the residual strength of 
the composite decreases due to impact damage, the stress wave factor 
decreases and the attenuation measured with the through transmission 
technique increases. Thus, a correlation between the SWF and attenuation 
was indicated also. 

In this study the output signal due to a single tone burst input 
will be analyzed theoretically and experimentally for the case where 
transmitting and receiving transducers are both coupled to the same face 
of the structure. 

If an ultrasonic disturbance is introduced into or is generated with- 
in an elastic isotropic structure, the parameters which govern its sub- 
sequent propagation and detection may be summarized as follows: 

1. Material Properties; Modulus of elasticity,\Poisson’ s 
ratio, density, shear and longitudinal attenuations, 
velocity, microstructural defect state,...; 

2. Geometrical properties; Overall geometry, discontinuities, 
macro structural defect state,...; 
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3. Measurement conditions: Location and size of transducers, 

sensitivity and frequency response of transducers, couplant, 
frequency and amplitude characteristics of input,... 

A more extensive list of parameters, including those which must be con- 
sidered for fiber reinforced composite structures is given by Williams [8]. 
Various authors have discussed the effects of material anisotropy and 
macrostructural defect state [9-11] as well as experimental measurement 
conditions including type and thickness of couplant [12-14] and transducer- 
structure contact pressure requirements [1], 

Finally, while it is clear that the proposed analysis has significant 
value in ultrasonic testing, it is important to reemphasize that such an 
analysis is equally valuable in acoustic emission. For example, since 
both longitudinal and shear waves may be generated by AE events and since 
both types of waves must propagate through the structure before being 
detected, it is clear that material properties, geometrical properties 
and measurement conditions all play a role in what is ultimately detected 
as the AE signal. Further, since the reflection at boundaries of one 
type of wave produces the other type, this indicates various coupling 
effects between the material, geometric and measurement features on both 
ultrasonic and acoustic emission parameters. 
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THEORETICAL ANALYSIS 


A schematic of the system under investigation is shown in Fig. 1. 

A transmitting and a receiving transducer of radius a and spaced a 
distance £ apart are coupled to the same face of an isotropic elastic 
plate of thickness h and of infinite planar (x - y) extent. The input 
electrical voltage to the transmitting transducer is V£(t) and the out- 
put electrical voltage from the receiving transducer is Vo(t) where t 
represents time. The output voltage Vo(t) due to a steady-state 
sinusoidal input voltage Vi(t) are analyzed. 


Assumptions on the Transducers 


The longitudinal transducers are assumed to transform an electrical 
voltage into a uniform normal stress and vice versa. The approach which 
follows is similar to [15]. 

If an input voltage of amplitude V and frequency (a) is applied 
according to 


V^(t) =76^“"^ , (1) 

the normal stress that is introduced into the specimen at the trans- 

ducer-specimen interface by the transmitting transducer is 


a (t) = F, (w)V e 
zz 1 


-i(oot + 


( 2 ) 


where F]^(oj) is the transduction ratio for the transmitting transducer in 
transforming a voltage to a stress and is a phase angle. In eqns. (1) 
and (2) , the harmonic character of the signals are expressed in the 
complex notation where i = and only the real parts of these and sub- 

sequent equations should be considered. Thus, the amplitude T of the 
applied stress is defined as 


T = F^(w)V (3) 

Similarly, if a stress wave producing a normal stress component 
Ozz of amplitude T’ and frequency (jJ that impinges on the receiving trans- 
ducer is defined as 
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(4) 
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the output voltage from the receiving transducer is 


-i(o)t + <!>„) 

V (t) = F_(w) e 
o 2 


(5) 


where F 2 (oa) is the transduction ratio for the receiving transducer in 
transforming a normal stress to a voltage and (p2 is a phase angle. Thus, 
the amplitude V’ of the output electrical voltage is 

V’ = F^Cco) T’ (6) 


Stress Field Radiated by the Transmitting Transducer 


The stress field radiated into an isotropic elastic semi-infinite 
solid by the transmitting transducer is now considered. Internal 
reflections from the faces of the plate are considered in a subsequent 
subsection. 

Fig. 2 shows an isotropic elastic half space bounded by the (x - y) 
plane and extending into the positive z-direction. The semi-infinite 
medium is excited by a uniform harmonic normal stress over a circular 
region of radius a at its surface. 

Using the cylindrical coordinates (r,4),z) as shown in Fig. 2, 
the boundary conditions at z = 0 require that the shear stresses 
and G 2 ;(P vanish and that the only nonzero stress is the normal stress 


a 

zz 


_ -loot ^ ^ . 

T e for 0 £ r £ a 

0 for r > a 


(7) 


This is an extension of the so-called Lamb problem [16] and using 
Hankel transform techniques. Miller and Pursey [17] have derived the 
displacement field solution at any point in the medium. From the dis- 
placements, the normal stress at any point in the half space can be 
shown to be 
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where 



(9) 


where is the longitudinal wave speed, C2 is the shear wave speed, ^ 
is an integration variable introduced in the Hankel transform, Jq and Jx 
are Bessel functions of the first kind and of order zero and one, 
respectively . 

The asymptotic behavior of eqn. (8) for large R and small a can be 
obtained using procedures in [17] as 


a - iaT 
zz 


Ik^CR - 


D^(o),e) 


Cit) 


ik^CR - 


402(10,6) 


C2t) 


( 10 ) 


where 

2 2 2 

cos 0 (k - 2 sin 0) J. (k a sin 0) 

- - 

sin 6 [ (k^ - 2 sin ^0 ) ^ + 4 sin^ 0 cos 0 J ^ ' - sin^6 ] 

( 11 ) 

cos 6 (2 - sin^ 0) - k^sin^ 6 Jj^(k 2 asin 0 ) 

D2 = — 

sin 0 [k(2 sin^ 0 - 1)^ + 4 sin^ 0 cos 6 ■/ 1 ~ k^sin^ 6 ] 
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and 


k 



The first term in the square bracket in eqn. (10) propagates with 
the longitudinal (P) wave velocity c^. Thus, this normal stress component 
is associated with the so-called P-wave. The second term is the square 
bracket in eqn. (10) propagates with the shear (S) wave velocity co. 

Thus, this normal stress component is associated with the so-called S-wave. 
The factors D]^ and D 2 are functions of frequency co and angle 0 and are 
the directivity functions for the P and S waves, respectively. Fig. 3 
shows polar diagrams of the directivity function associated with the 
P-wave for various frequencies. 


Near Field and Far Field 


The asymptotic solution given by eqn. (10) is valid for large R or 
the far field. In order to determine the magnitude of R at which the 
asymptotic solution becomes valid, the truncated terms in the asymptotic 
evaluation of eqn. (8) should be investigated. However, by using an 
analogy with a similar acoustics problem, an estimate for the far-field 
distance can be obtained readily. 

The acoustics problem considered is a rigid piston vibrating in a 
rigid baffle. The radiation pressure along the piston axis in the 
acoustic medium is of interest and is given in [18]. The amplitude of 
the radiation pressure fluctuates near the piston and eventually 
approaches the asymptotic solution given in [19]. The near-field 
distance Zq within which rapid pressure fluctuations occur is [20] 


z 

o 



( 12 ) 


where X is the wavelength. Then the far field distance zp where the 
asymptotic solution applies is given as [2] 


Zp = 3z^ (13) 

Combining eqns . (12) and (13), the P-wave component of 0^2 given 
in eqn. (10) is valid if 


R > R^ =3 


2 2 
4a - 


4A 


1 


(14) 
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( 15 ) 


where is the far-field distance and is the longitudinal wave- 

length and is given by 


X 


1 


2ttc 


1 


03 


Wave Reflection at a Plane Boundary 


Waves in the plate specimen experience multiple reflections between 
the faces of the plate. An obliquely incident plane P-wave is reflected 
from the boundary as a P-wave with an angle of reflection equal to the 
angle of incidence, say and with a S-wave with an angle of reflection, 

say 02 , where 0 ]^ and 02 are related by Snell’s law as [ 9 ] 


^1 

sin 0, = — sin 0„ ( 16 ) 

1 C 2 2 

The magnitudes of the reflected P and S waves are defined by the reflection 
coefficients Qpp and Qpg, respectively, and are given by [ 9 ] 

2 2 

sin 20^ sin 20„ - k cos 20„ 

Q = ^ ( 17 ) 

^P? 2 2 

sin 20^ sin 262 + k cos 2©2 

and 

2 sin 20- cos 20„ 

Qps i ( 18 ) 

sin 20j^ sin 202 + k cos 202 

Spherical waves are generated in the plate specimen shown in Fig. 1 . 
However, incident and reflected spherical waves are nearly plane near the 
boundary of reflection if the wavelengths are much smaller than the 
distance travelled by the wave from the wave source to the boundary [21]. 
Then, the plane wave reflection solution given in eqns. ( 16 ) through ( 18 ) 
applies . 


- 7 - 



Steady-State Output Voltage Amplitude due to Multiple Wave Reflec tions 
in a Plate 


Fig, 4 illustrates a plate of thickness h with a transmitting trans- 
ducer located at point 0 and a receiving transducer located on the same 
face at point M, separated by a distance £. A path in the plate is also 
shown where the wave is reflected n times from the bottom face of the 
plate. Only P~waves are considered in what follows. The total number of 
reflections experienced by the wave in travelling from the transmitting 
transducer to the receiving transducer is (2n - 1). With respect to z 
the angle of incidence of the wave at either face of the plate is 0]^ and 
the total distance travelled by the wave is From the geometry in 

Fig. 4, 


0 


1 


-1 

tan 


^2nh^ 


(19) 


and 


R 

n 


I 

sin 0j^ 


( 20 ) 


The time delay tj^ for the wave to reach the receiver after n reflections 
from the bottom face is 


t 

n 



( 21 ) 


If there were no bottom boundary and the waves were propagating in 
an infinite half space, the waves would travel to point M' in time t^- 
Then the amplitude of the hypothetical normal stress at point M’ is 
T^* and is defined per eqn. (10) as 


T 


M 


f 


aT 

n 


( 22 ) 


However, with the bottom boundary present, the wave is reflected a total 
of (2n - 1) times. Thus the amplitude of the normal stress at point 
M is Tj 4 and is obtained by modifying eqn. (22) as 
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( 23 ) 


T = 
M 


R 


where Qpp is given in eqn. (17). The amplitude of the output voltage 
from the receiving transducer is Vo and can be obtained by combining eqns. 
(6) and (23) as 


V 

n 


a 

n 



1 


(24) 


Substituting eqn, (3) into eqn. (24) gives 


V 

n 


V a 



R 

n 


(25) 


Introducing the longitudinal wave attenuation constant a of the medium 
and a possible signal amplification factor K, eqn. (25) can be rewritten 
as 


V 

n 


K V a 



1 -aR^ 

e II 

R 

n 


(26) 


Eqn. (26) gives the output voltage amplitude from the receiving transducer 
due to an input voltage amplitude V at the transmitting transducer when 
the wave path has included n reflections from the bottom face of the 
plate specimen. Only P-waves have been considered and eqn. (26) is valid 
if the plate thickness h is larger than the far-field distance Rp^. 



EXPERIMENTAL EQUIPMENT AND PROCEDURES 


A schematic of the experimental system is shown in Fig, 5, The 
system consisted of a pulsed oscillator (Arenburg model PG-652C) for 
generating the sinusoidal waves; a low frequency inductor (Arenburg 
model LFT-500) ; a low frequency amplifier (Arenburg model LFA-550) ; broad- 
band (0.1 to 3,0 MHz) transmitting and receiving transducers (Acoustic 
Emission Technology CAET) FC-500) having an approximately flat sensitivity 
of -85 dB (relative to IV/yBar) ; a transducer specimen interface couplant 
(AET SC-6); and an oscilloscope (Tektronix model 455). An attenuator, 
set at 10 dB, reduced the input signal to 100 V (peak-to-peak) into the 
transmitting transducer. 

The specimen was an aluminum (6061-T6) plate 10 cm in thickness and 
25 cm on a side. The plate specimen was supported by a 2.5 cm wide 
rubber mount on the bottom face. The diameter of the wearplate of the 
FC-500 transducer was 2.25 cm and the clamping pressure on the transducer 
exceeded the saturation pressure which is defined in [1] as the minimum 
transducer-specimen interface pressure which results in the maximum 
amplitude of the output signal, all other parameters being held constant. 

A single 100 V peak-to-peak tone burst (finite duration sinusoidal 
signal) was input into the transmitting transducer. The output from the 
receiving transducer was recorded at various spacings between the trans- 
ducers and tone burst input frequencies. The duration of the tone burst 
input was chosen so that the output signals due to multiple reflections 
did not overlap. The input frequencies were 0.75, 1.50 and 2.25 MHz. 

The transducer spacings (il) between the transducer axes were 5.5 and 
10,0 cm. 

The aluminum had a density of 2.7 g/cm^; longitudinal wave velocity 
of 6.32 X 105 cm/sec [2]; shear wave velocity of 3.13 x 10^ cm/sec [2]; 
surface or Rayleigh wave velocity of 2.93 x 10^ cm/sec [9]; and longitudinal 
wave attenuation of 0.015 neper/cm [22]. Note that it was assumed that 
the longitudinal wave attenuation was constant for the frequencies 
considered. 
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RESULTS AND DISCUSSIONS 


Typical input and output signal records are shown in Fig. 6. The 
input signal frequency was 1.50 MHz and the transducer separation was 
5.5 cm. The time delays for multiple reflections are indicated both 
above and below the output .record . The notations P, S and R represent 
the longitudinal, shear and Rayleigh waves, respectively. Also the 
notations P-P and P-S represent a P-wave reflected from the bottom 
face producing a P and a S wave, respectively. 

Fig. 6 shows that the output signal is dominated by reflections of 
P-waves only. (Refer to the time delays above the output record.) 

Because the specimen is finite in its planar dimensions, reflections 
from the vertical faces began to appear after the 6th reflection from the 
bottom face. Fig. 6 also shows that for this special case, the 2nd 
reflection from the bottom face produces the output wave packet having 
the maximum amplitude. 

Table 1 shows the calculated and measured time delays t^ for the 
P-waves to reach the receiving transducer after multiple reflections 
from the bottom face. The calculated values are based on eqn. (21). 

The measured values are accurate to within 1 ]Js. Table 1 shows that 
there is excellent agreement between the calculated and measured time 
delay values. 

Theoretical results based on eqn. (26) are plotted with the 
experimental results in Figs. 7 through 10. The vertical axis represents 
the amplitude of the output signal of each reflection from the bottom 
face. The horizontal axis is the time scale indicating the time delay 
for each reflection to travel from the transmitting transducer to the 
receiving transducer. Only P-waves are considered. Because the factor 
FiF 2 in eqn. (26) is not known, the theoretical result is arbitrarily 
specified to match the experimental result at the 4th reflection. (In 
fact, this may be considered a technique for estimating F 1 F 2 .) Fig. 7 
is obtained at 0.75 MHz with a transducer spacing of 5.5 cm; Fig. 8 is 
obtained at 1.50 MHz with a transducer spacing of 5.5 cm; Fig. 9 is 
obtained at 1.50 MHz with a transducer spacing of 10.0 cm and Fig. 10 
is obtained at 2.25 MHz with a transducer spacing of 5.5 cm. 

There is good agreement between the theoretical and experimental 
data in Figs. 7 to 10. In every case, the maximum discrepancy between 
the theoretical and experimental data occurs for the first reflection 
from the bottom face. The probable reasons for this discrepancy are 
one or more of the following: 

1. The theoretical analysis assumes that the receiving trans- 
ducer produces an output voltage that is proportional to 
the normal stress applied to its face. However, only the 
normal stress at the center of the receiving transducer 
is used in the study. The normal stress varies over the 
transducer face depending on the angle of incidence, the 
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wavelength and the transducer diameter. Given the same 
wavelength and the same transducer diameter within a test, 
the maximum error will occur when the angle of incidence 
is largest; and this happens for the first reflection. 

2. The theoretical study considers a steady-state input with 
a single frequency. However, the tone burst input used in 
the experiments has a finite duration. Thus, the input 
frequency spectrum has a finite bandwidth. Yet, only the 
center frequency of the spectrum is considered in the 
theoretical calculations. From Fig. 3, it is observed 
that at large angles of incidence, the directivity function 
gives decreasing values of amplitude for increasing 
frequencies. Thus for the first reflection where the angle 
of incidence is largest, the error is expected to be largest 
compared with other reflected waves because the higher 
frequency content of the spectrum has been drastically re- 
duced by the directivity function. This results in the 
greater distortion of the signal recorded by the receiver 
for the first reflection as shown in Fig. 6. Furthermore, 
these arguments suggest that the theoretical prediction 

of the amplitude of the first reflected wave should exceed 
the experimentally measured amplitude. This is consistent 
with Figs. 7 through 10. 

3. The theoretical solution in eqn. (26) is based on an 
asymptotic approximation for distances larger than the far- 
field distance . The far-field distance for 0.75, 1.50 

and 2.25 MHz are calculated using eqn. (14) to be 3.89, 

8.73 and 13.35 cm, respectively. The thickness of the plate 
specimen is 10.0 cm. Because the distance travelled by the 
first reflected wave is less compared to other reflected 
waves, the error in the asymptotic theoretical solution 
is expected to be larger for this first reflected wave. 
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CONCLUSIONS 


Ultrasonic input-output characteristics for an isotropic elastic 
plate with transmitting and receiving longitudinal transducers coupled 
to the same face have been analyzed theoretically and experimentally. 

In the theoretical analysis, it was assumed that the transmitting 
transducer applied a uniform normal stress to the specimen surface when 
the transducer was excited by an electrical voltage. The normal stress 
field in the plate was calculated via a boundary value problem for an 
isotropic elastic semi- inf inite medium subjected to a uniform normal 
stress applied to a circular region at the surface and varying 
harmonically in time. An asymptotic solution for the normal stress 
field was obtained and the frequency and geometric dependent directivity 
functions for the radiated stress field were found. The requirement 
for the asymptotic solution to be valid was also discussed. Then, the 
radiated stress waves were confined to the plate specimen by considering 
wave reflections at the top and bottom faces. Reflection coefficients 
were used to describe the wave reflections behavior at the faces. By 
introducing material attenuation and by assuming that the receiving 
transducer transforms the incident normal stress into an electrical 
voltage, the output voltage amplitude of the receiver was derived. 
Particular attenuation was given to the P-waves only. 

There was good agreement between the theoretical and experimental 
results for the output voltage amplitudes due to multiple reflections 
of longitudinal waves. Similarly, excellent agreement between the 
theoretical and experimental results for the times-of-arrival of the 
individual wave packets at the receiving transducer were obtained. 

This study provides a major step forward in the quantitative 
understanding of ultrasonic nondestructive evaluation (NDE) parameters 
such as the stress wave factor (SWP) . It also provides the potential 
for assisting in the development of more efficient and more revealing 
NDE schemes utilizing the SWF. Acoustic emission, ultrasonic attenuation 
measurements as well as transducer calibration techniques may also 
benefit from this study. Finally, this work should be extended to 
the interrogation of composite materials. 
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TABLE 1 Calculated and Measured Time Delays for Longitudinal 

Waves to Reach the Receiving Transducer after Multiple 
Reflections 


Number of 
Reflections from 
Bottom Face 
n 

Time Delay t , (ys) 
n 

for £ = 5 

5 cm 

for £ = 10 

0 cm 

Calculated 

Measured 

Calculated 

Measured 

1 

32.8 

32 

35.4 

35 

2 

63.9 

63 

65.2 

65 

3 

95.3 

95 

96.3 

96 

4 

126.9 

126 

127.6 

127 

5 

158,5 

157 

159.0 

158 

6 

190.1 

188 

190.5 

188 

. 
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Input Signals: 




Fig. 1 Schematic of ultrasonic (stress wave factor) configuration showing typical input signals. 




Fig. 2 Half-space problem in evaluating the radiated stress field due 
to the transmitting transducer in the cylindrical coordinate 
system (r,4),z). 


- 18 - 



0.75 MHz 


.50MHz 


2.25MHz 


Polar diagrams of frequency and angle dependent directivity function Di for the normal stress 
due to longitudinal waves with transducer radius a of 1.125 cm, longitudinal wave speed cj^ 
of 6,32 X 10^ cm/sec, and shear wave speed C2 of 3.13 x 105 cm/ sec. 





Pulsed Oscillator 


Oscil loscope 



Fig. 5 Schematic of experimental system for measuring input-output for transmitting and receiving 
transducers coupled to the same face of a plate specimen. 





Fig. 6 Input signal (a) to the transmitting transducer and the 

corresponding output signal (b) from the receiving transducer 
at 1.50 MHz and a separation between transducers of 5.5 cm. 
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